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The framework of Large Extra Dimensions provides a way to explain why gravity 
is weaker compared to the other forces in nature. A consequence of this model is the 
possible production of D-dimensional Black Holes in high energy p-p collisions at the 
Large Hadron Collider. The present work uses the CATFISH Black Hole generator to 
study quantitatively how these events could be observed in the hadronic channel at 
mid-rapidity using a particle tracking detector. 
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1. Introduction 

With the Large Hadron Collider (LHC) opening in late 2008, physics at the 14 TeV 
center of mass scale will be probed in p-p collisions. Although searching for the 
Higgs particle will be one of the main research focuses, other interesting physics, 
such as the possible formation of D-dimensional Black Holes (BH) in these collisions 
will be another direction of intense research. This phenomenon is predicted by the 
framework of Large Extra Dimensions (LED) in the Arkani-Hamed, Dimopolus and 
Dvali (ADD) model^^l This model is a proposed solution to the Hierarchy Problem, 
the question of why gravity is much weaker than the other forces in nature. The 
scale for gravity is characterized by the Planck mass, Mp ss 10 19 GeV, whereas the 
strong and the weak forces have a scale on the order of 1 GeV/fm and 100 GeV, 
respectively. 

The approach to solve the problem in the ADD model uses the following four 
assumptions: 1) the hierarchy of the forces in nature only exists in 3+1 dimensions, 
but not in higher dimensional space-time with D=3+l+n dimensions, 2) only grav- 
ity is allowed to propagate in the n extra dimensions through gravitons, 3) the 
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LED are "compact" or finite, but too small to be detected until (possibly) now, 4) 
Standard Model particles only propagate in 3+1 dimensional space-time, which is 
embedded in 3 + 1 + n higher dimensional space-time. Gravity seems to be weaker 
in 3+1 dimensions because it is diluted in 3+1+n higher dimensional space. The 
Planck mass in a theory with n LED is given by 



where Mpo = y/h/G = 1.22 x 10 19 GeV is the Planck mass in D=3+l dimensions 
and r c is the compactification length of the extra dimensions. The compactification 
length is a free parameter of the theory, and it can be seen El that for a given Planck 
mass, it decreases as the number of extra dimensions increases, for example for Mp 
= 1 TcV r c is close to 1 fm for n=7. For fixed r c the Planck scale is lowered as the 
number of LED increases and the hierarchy of the forces is largely suppressed 

One of the consequences of this model is the formation of D-dimensional BHs 
smaller than the size of the extra-dimensions and centered on the brane. The grav- 
itational radius of the Z?-dimensional BHs is up to 10 32 times larger than that of a 
usual BH in 3+1 dimensions with the same mass. This considerably increases the 
possibility of creating such an object at the LHC. 

The main purpose of this work is use the CATFISH SI Monte Carlo (MC) code 
to study the possible signatures of BHs formed in p-p collisions, which could be 
observed in the hadronic channel at mid-rapidity using a charged particle tracking 
detector. CATFISH is a MC code based on the ADD model and is used to generate 
BH events at the LHC with center of mass energy of 14 TeV. It is similar to other 
MC generators like TRUENOIR S and CHARYBDIS but it is more up to date 
since it uses the most recent theoretical results of BH formation and evolution. 
The updates include inelasticity effects during the BH formation phase, exact field 
emissivities, corrections to Hawking semi classical evaporation phase, BH recoil on 
the brane, and additional final BH decay modes (including remnants). CATFISH 
links to the PYTHIA ^ MC code to simulate the evolution of the decay products 
of the BH into Standard Model particles. The flexibility of CATFISH allows the 
study of the different theoretical models of BH formation. A comparison between 
the features of CATFISH and CHARYBDIS has been carried out before ^ and a 
similar study of the observables of BH at the LHC has also been done using the 
CHARYBDIS Code^. 

CATFISH includes the following three models to calculate the BH evolution: 

• No Gravitational Loss model (NGL): This model works in the semi 
classical limit of BH formation. The hoop conjectured^ is used to estimate 
the possibility of BH formation in a particle collision. It states that an 
apparent horizon is formed in D=A dimensions if a mass M is compacted 
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in a region with a circumference C such that 

where r^M) is the Schwarzschild radius for the mass M given by 

fh{M) = — -g- . (3) 



{D-2)n D _ 2/ 

Here Qd-2 is the volume of the (D — 2)-sphere and Gu is the gravitational 
constant. The impact parameter b of two colliding partons has to be smaller 
than rh{M) to produce a BH. The cross section is given approximately by 
the geometrical Black Disk (BD) cross section 

OBD = irR 2 {s,n)<d[R{s,n) - b] (4) 

where R is the horizon radius and depends on the center of mass energy y/s, 
and the number of extra dimensions n. The BH mass is equal to the center 
of mass energy of the partons forming the BH. This is the same model used 
in TRUENOIR and CHARYBDI S. 

Yoshino-Nambu model (YN) This model uses the Trapped Sur- 

face approach, which gives a bound on the inelasticity of a collision by 

1 1 Ol 

modeling two incoming partons as Aichelburg-Sexl shock waves^. The ap- 
parent horizon is found in the union of the two shock waves. The condition 
for BH formation is better described in higher dimensionalities by the vol- 
ume conjecture than by the hoop conjecture, and is given by 

l/OD-3) 



Vn- 3 



< 1, (5) 



where fi.D-3 is the volume of the (D — 3)-sphere and Vd-3 is the charac- 
teristic (D — 3)-dimensional volume of the system. The volume conjecture 
reduces to the hoop conjecture in D = 4 dimensions. The cross section is 
calculated in this model as 

PBHproduct ion = F \D)tTt1{s , Tl) , (6) 

where F(D) is a numerical factor close to unity. The BH mass is less than 
the center of mass energy of the partons forming the BH due to emission 
of gravitons and it depends on the impact parameter. 



• Yoshino-Rychkov mo del (YR) 031 This model is an improved version 
of the YN model in that the apparent horizon is constructed from a slice of 
the future light cone in the shock collision plane. The slice is to the future 
of that used in the YN model. The condition for BH formation is also given 
by the volume conjecture and the cross section calculation is similar to the 
calculation for Equation [6l The BH mass is also reduced in comparisson to 
the BH mass in the NGL model. 
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The event simulation in CATFISH occurs in three steps: 

(i) The initial BH mass is sampled from the differential cross section. 

(ii) The BH is decayed through the Hawking mechanism and final hard events. 

(iii) The unstable quanta emitted are hadronized or decayed instantaneously by 
PYTHIA, except top quarks, which are decayed as t — > bW first. 

The results of our calculations are presented in the following section. 



2. Results and Calculations 

The signatures of BH production are studied in this section, first for completely 
decaying BHs and then for BH remnants. The results of our calculations are shown 
in Figures 1-12. 



2.1. BH signal using the hadronic channel 

The hadronic channel is considered as a possible method to detect BHs. To be 
conservative, the first-year luminosity, the transverse momentum resolution and the 
rapidity are assumed to be L = 10 31 cm _2 s _1 , 0.1 GeV/c < pr < 300 GeV/c and 
-1 < y < 1 respectively (e.g. these characteristics are similar to those of the central 
tracking detectors of the ALICE experiment ^ a t CERN) . 

The transverse momentum distribution — 4^- ^^BH f or BH events flattens more 

Pt dp T 

than in the QCD background events at pr > 200 GeV/c for Mp=l TeV and at 
higher pp for larger values of Mp. This would make the detection of BH events 
difficult for tracking detectors since the momentum resolution of these detectors is 
poor for large px- To overcome this problem another observable can be used instead, 
namely the sum of the transverse momentum of all charged hadrons in each event, 
calculated using 



P T = J2 PTt , (7) 

i 

where i runs over all charged hadrons in one event. CATFISH is used to generate 
BH events and for each the Pt is calculated and histogrammed in bins of size of 0.1 
TeV/c. The total number of counts per bin for an arbitrary period of time N can 
be calculated using the following relation: 



Nhi n <jLt . . 

N= _Jnn ^ (8) 

J^Tot 

where Nbin is the number of counts in a particular bin, Npot is the total number of 
events generated, and a is the cross section. For the each simulation in this work, 
Nrot = 10 5 events were generated and a period of time of four months was chosen, 
so that t w 10 7 s. 
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In Figure \\\N vs. Pt is plotted for the NGL model with different values of the 
Planck mass, together with a normal QCD background run from PYTHIA. This run 
is a mixture of six different runs changing the hardness of the 2 — > 2 parton collision. 
It is seen that at low Pt the normal QCD events dominate, but above 0.5 TeV the 
Mp = 1 TeV curve takes over. This is a possible method for detecting the BHs: i.e, 
noticing the discontinuous change in the slope from normal QCD background 
events to BH events. For higher values of the Planck mass the curves intersect the 
QCD background at higher Pt values. The BH production becomes more suppressed 
but it can be seen that even at a Mp of 5 TeV it would still be possible to detect 
a signal under these running conditions. Furthermore, by obtaining experimentally 
the Pt value at which the change in slope happens, the actual value of the Planck 
mass could be determined. 

Figure [2] plots the same quantities as in Figure [1] except that various cuts are 
applied to simulate the acceptance of a realistic tracking detector. The cuts used are 
on rapidity (—1 < y < 1), transverse momentum ( 0.1 GeV/c < pt < 300 GeV/c) 
and multiplicity in the rapidity cut (m >65). As seen the BH signal is not degraded 
by these cuts. Its effect is particularly visible in lowest Pt bins where the number 
of counts is reduced. These cuts will be applied in all subsequent plots. 

Figure [3] plots the combined signal of the QCD background and the BH con- 
tribution, as one would obtain from the experiment. It can be seen that at Pt ~ 
0.5 TeV/c, the slope of the graph changes abruptly indicating the transition from 
normal QCD events to BH events. 

In Figure [3] the effect of varying the number of extra dimensions is studied. As 
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Fig. 1. Comparison of the Pt for different values of Mp with the QCD background. Model: NGL 
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can be seen there is no substantial difference in our res ults if five, six or seven extra 
dimensions are used in the simulations. String Theory ^2 favors n=7 , so we have 
chosen to use this value for all of the results presented in the rest of the study. 
In Figure [5] the Pt is plotted for the three theoretical models of BH formation: 
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Fig. 2. Comparison of the with applied cuts for different values of Mp and the QCD back- 
ground. Model: NGL. 
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Fig. 3. Combined signal of the QCD background and the BH signal for a Mp = 1 TeV. Model: 
NGL. 
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NGL, YN and YR. The NGL model has the largest cross section of the three models 
since it uses the semi classical approximation to BH formation. The YN and YR 
models take into account the gravitational energy loss at formation and this de- 
creases the cross section. It can be seen that only the NGL curve is above the QCD 
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Fig. 4. Comparison of the Pt for Mp of 1 TeV for different number of extra dimensions and the 
QCD background. 
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Fig. 5. Comparison of the three different BH formation models for Mp = 1 TeV: No Gravitational 
Loss (NGL), Yoshino-Nambu (YN) and Yoshino-Rychkov(YR). 
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background, meaning that if BHs are formed according to the YN or YR models 
it would not be possible to detect any using the Pt method. From this figure we 
assert that the NGL model is an upper bound to the formation of BHs, and the YN 
and YR are the lower bounds for the known models. 

In the p-p collisions BHs are formed when a given minimum mass is reached. 
This can be varied in CATFISH using the parameter X m i n > 1 by 

M mm ( formation) = X min x M P , 

Subsequently BHs evaporate through the Hawking mechanism until a minimum 
mass is reached, 

M min (evaporation) = Q m in X Mp 

with Qmin > 1. It is assumed in our simulations that the minimum BH mass at 
formation and at evaporation are the same, but in general these two values are 
independent. Figure [5] shows results for how the minimum BH mass at formation 
and at evaporation affect the hadronic signal. The Pt is plotted for a Mp of 1 TeV 
and two values of X min — Q m i n using the NGL model. It is observed that at low 
Pt the BHs with lower initial mass have a higher signal. This is because the decay 
product of the BHs with larger masses would have a higher Pt than that of BHs 
with lower masses. For values above 1.5 TeV/c the Pt is about the same for both, 
showing that BH production is insensitive to this parameter for higher Pt- 

So far it has been shown that the possibilities of detecting a BH using the 
sum of the transverse momentum depend strongly on the model of BH formation 
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Fig. 6. Comparison of two values of the BH mass at formation (Xmin) and BH mass at evapo- 
ration (Qmin) and the QCD background. Model: NGL with M P = 1 TeV . 
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and evaporation. The NGL model provides optimistic chances of producing BHs, 
whereas the YN and YR models provide a more pessimistic expectation, because 
the BH events would not be seen above ordinary QCD background. Note that this 
does not imply that BHs cannot be detected in particle trackers at mid-rapidity. In 
the next section we explore such a possibility by studying the detection of charged 
remnants, if they are formed. 

2.2. Signal from BH Remnants 

The possibility of a BH not decaying completely to Standard Model particles after 
the Hawking evaporation phase and leaving a remnant has also been discussed in 
the literature 1 ^. The BH remnant would be a very massive and possibly charged 

particle and these facts could be used to detect them in mid-rapidity particle 
trackers. The remnants are produced in CATFISH when the evaporation phase is 
over and the mass has reached the value Mbhr = Qmin x Mp, which is not further 
decayed into Standard Model particles. 

Using the technique of summing the transverse momentum it is seen in Figures 
17181 and [5] that having a BH remnant significatively reduces the total number of 
counts per bin. Even in the most optimistic scenario, the NGL model, Figure [7] 
shows that the BH remnant events are below the QCD background. Therefore the 
possibility of detecting the BHs via the Pt method are very small because the QCD 
background dominates. 

However, in light of having a BH remnant, anot her technique could be used to 
detect them. BH remnants may have a net charge and that would allow one 
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Fig. 7. Comparison of the Pt in the NGL model with no remnant and with BH remnants of 1 
and 2 TeV 
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to detect these massive particles using the time-of-flight (TOF) method in particle 
trackers. By knowing the momentum of the BH remnant and the TOF, its mass 
can be reconstructed using the following relation: 
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Fig. 8. Comparison of the Pt in the YN model with no remnant and with BH remnants of 1 and 
2 TeV. 
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Fig. 9. Comparison of the Pt in the YR model with no remnant and with BH remnants of 1 and 
2 TeV. 
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TP 



(9) 



where p' is the momentum of the remnant. (3' and 7' can be calculated respectively 
from 



3' = -<y'= 1 



(10) 



We assume a time resolution of 50 ps, a momentum resolution Ap/p I0%p and 
a flight path of x = 3 m. In Figure [TOl we plot the counts per 0.05 TeV/c 2 bin of 
the actual BH remnant mass obtained from the simulation and the reconstructed 
mass from the TOF and momentum. It can be seen that both graphs are almost the 
same for this bin size, thus demonstrating that the reconstructed mass is correct. 
The BH remnant mass is observed to have a distribution and not the exact value 
that was used in the simulation Q m in x Mp. As will be seen in Figure fT3l the mass 
distribution becomes broader as the value of the minimum BH mass at evaporation 
gets closer to the Planck mass. 

In Figure [IT] we plot the difference between the actual mass of the BH and the 
reconstructed mass using a 0.2 GeV bin size. The mass difference has an average 
of GeV and a standard deviation of 0.76 GeV for this particular case, showing 
the high accuracy which is possible in determining the BH mass using a charged 
particle tracker plus TOF with reasonable performance characteristics. 
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Fig. 10. Comparison of BH remnant mass with the reconstructed mass from the TOF for Mp = 
1 TeV, Q min = 1 
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CATFISH does not assign a charge to the BH remnants and assumes they are all 
neutral, t her efore in our simulations we assign a charge to each remnant, following 
reference The charge distribution among the BH remnants is plotted in Figure 
[T2l It is approximately Gaussian with a standard deviation of 1.47. The plot shows 
that about 75 percent of the BH remnants will have a net charge and could be 
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detected using the tracking detector plus the TOF technique. 

As was mentioned earlier in the paper, the BH remnant mass has a distribution 
which is broader when it is closer to the Planck scale and it gets narrower when it is 
a few times larger. Figure [T3l shows a plot of how the BH remnant mass distribution 
changes depending on how high the minimum mass is above the Planck scale. One 
of the strong points of this method of detecting BHs is that the signal should be very 
clean in these large mass regions since no QCD process can generate such masses 
so that even small signals should be recognizable. 

3. Summary 

We have studied some of the possible methods which could be used to detect BHs 
at the LHC, if they are formed, by looking in the hadronic channel at mid-rapidity. 
Finding the sum of the total transverse momentum in charged hadrons is one such 
method, because it is expected that the decay products of BH events have larger 
Pt bin counts than normal QCD events at Pt > 0.5 TeV in the NGL model. The 
BH signature here is the transition from background QCD events to BH events. 
This model represents the upper bound for BH formation because there is no grav- 
itational energy loss during formation. The other two models studied, the YN and 
YR models, take into account the energy loss at formation, and thus provide the 
lower bounds for the known models of BH formation. In this case lower counts 
per Pt bin are expected and the hadronic signal will not be recognizable from the 
QCD background events. If the Pt method does not work, another way of detecting 
BHs is possible if they do not decay entirely and, instead, leave a remnant. Using 
tracking detectors plus the TOF, the mass of a BH remnant can be reconstructed. 
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Comparison of different BH remnant mass distributions for Mp = 1 TeV. 
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By looking at the spectrum of masses from many events it could be easily recog- 
nized that there is a very massive particle left over. This would be the signature 
for BH events with remnants happening at the LHC. There are other possible BH 
signatures which h ave not been considered here and have been studied elsewhere 
| 4 | 19 | 20 | 21 | 22 | 23 | 2 4] Examples include the missing energy, missing transverse energy 
and hadron energy signatures, as well as signatures due to suppression of back-to- 
back-correlated di-jets and di-lepton production with large transverse momentum 



Acknowledgments 

The authors wish to thank Professor Marco Cavaglia for his kind assistance in help- 
ing us run his code CATFISH. We also wish to thank Eric Anderson for corrections 
to the text and D.G Humanic for helpful suggestions. This work was supported by 
the U.S. National Science Foundation under grant PHY - 0653432. 



References 

1. N. Arkani-Hamed, S. Dimopoulos and G. Dvali, Phys.Lett.B429, 263(1998). 

2. N. Arkani-Hamed, S. Dimopoulos and G. Dvali, Phys.Rev.D59,086004(1999). 

3. Thomas J. Humanic, Extra-dimensional physics with p+p in the ALICE Experiment, 
ALICE Internal Note, ALICE-INT-2005-017 version 1.0 (2005) 

4. M. Cavaglia, R Godang, L. Cremaldi and D. Summers, "CATF ISH: A Monte Carlo 
simulator for black holes at the LHC" , arxiv: | |hep-ph/0609001f 2 (2007). 

5. S. Dimopoulos and G. Landsberg, in: Proc. of the APS/DPF/DPB Summer Study on 
the Future of Particle Physics (Snowmass 2001) ed. N. Graf, eConf C010630 (2001) 
P321. 

6. C. M. Harris, P. Richardson, and B. R. Webber, [ar"xiv:hep-ph/0409309| ^l (2004). 

7. T. Sjostrand, L. Lonnblad adn S. Mrenna, PYTHIA 6.2 Physics and Manual, 
|arXiv:hep-ph/0108264| (2001). 

8. Douglas M Gingrich, Comparison of Black Hole Generators for the LHC, 



arXiv:hep-ph/0610219vl (2006). 



9. Thomas J. Humanic et al., Signatures for Black Hole Production from hadronic ob- 
servables at the Large Hadron Collider, [arXiv:hep-ph7u 607097yl (2006). 

10. K.S. Thorne is Magic without Magic: John Archbald Wheeler, edited by J. Klauder 
(Freeman, San Francisco, 1972). 

11. H. Yoshino and Y. Nambu, Phys.Rev. D 67 024009 (2003). 

12. M. B. Voloshin, Phys. Lett. B 518 (2001) 137; M. B. Voloshin, Phys. Lett. B 524 
(2002) 376 [Erratum-ibid. B 605 (2005) 426]; T. G. Rizzo, JHEP 0202 (2002) 011; V. 
S. Rychkov, Phys. Rev. D 70 (2004) 044003. 

13. O. I. Vasilenko, "Trap surface formation in high-energy black holes collision," 
|arXiv:hep-th/0305067l 

14. H Yoshino and V. S. Rychkov, Phys. Rev. D 71 (2005) 10402, 

S. B. Giddings and V. S. Rychkov, Phys. Rev. D 70 (2004) 104026. 

15. ALICE Collaboration, J. Phys. G: Nucl.Part.Phys. 30 (2004), 1517; ibid, 32 (2006) 1295. 

16. E. Witten, arXiv:hep-th/212247vl (2002). 

17. Benjamin Koch et al., Black Hole remnants at the LCH, arXiv:hep-ph/0507138 v2 
(2005). 



November 15, 2008 23:27 WSPC/INSTRUCTION FILE BHPaperRevi- 
sionl 



Quantitative Calculations for mini-Black Hole production at the Large Hadron Collider 15 

18. S.Hossenfelder , B.Koch, M. Bleicher, Trapping Black Hole Remnants 
|arXiv:hep-ph/0507140| (2006) . 

19. B. Koch, M. Bleicher and S. Hossenfelder, JHEP 0510 (2005) 053. 

20. H. Stocker, "Stable TeV - black hole remnants at the LHC: Discovery through dijet 
suppression, mono-jet emission and a supersonic boom in the quark-gluon plasma," 
|arXiv:hep-ph/ 0605062 

21. G. L. Alberghi et al., "Probing quantum gravity effects in black holes at LHC," 
|arXiv:hep-ph/060l"243| 

22. J. Tanaka, T. Yamamura, S. Asai and J. Kanzaki, Eur. Phys. J. C 41 (2005) 19. 

23. L. Lonnblad, M. Sjodahl and T. Akesson, JHEP 0509 (2005) 019. 

24. Arunava Roy, Marco Cavaglia, Discriminating Supersymmetry and Black Holes at the 
Large Hadron Collider arXiv:hep-ph/0801.3281(2008) 



